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SPACE -WEIGHTED COMMUNICATION PATH EST IMATIgltL^ 

The present inventxon relates to a method of 
estimating a communication path formed by a plurality of 
channels. Thus it relates to a technique referred to as 
5 reception diversity, whereby a receiver has a plurality 
of antennas each associated with a different 
communication channel. In other words, the invention 
proposes a method of estimating the impulse responses of 
the transmission channels. 

10 In a communications system, especially a radio 

communications system, the receiver receives for each 
communication channel a signal transmitted by a 
transmitter. The transmitted signal is subject to 
amplitude and phase fluctuations in the communication 

15 channel with the result that the signal received by the 
receiver is not identical to the transmitted signal. 
Signal fluctuations are essentially due to what the 
skilled person refers to as intersymbol interference. 
This interference can result from the modulation law 

2 0 employed for transmission and is also caused by multipath 
propagation in the channel. 

It is found that the received signal is generally 
the result of a large number of reflections in the 
channel. The various paths taken by the transmitted 

2 5 signal cause various delays at the receiver. The impulse 
response of the channel represents all such fluctuations 
affecting the transmitted signal. It is therefore the 
fundamental characteristic representative of transmission 
between the transmitter and the receiver. 

30 The impulse response of the channel is used in 

particular by an equalizer whose precise function is to 
correct intersymbol interference in the receiver. A 
standard method of estimating the impulse response 
consists in placing a training sequence made up of known 

35 symbols in the transmitted signal. 

The sequence is chosen as a function of the 
modulation law and the dispersion of the channel. In the 
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present context, "dispersion" is to be understood as 
meaning the delay affecting a transmitted symbol taking 
the longest path of the channel relative to the same 
symbol taking the shortest path. The dispersion is 
5 routinely expressed as a multiple of the time between two 
successive transmitted symbols, i.e. a number of "symbol 
periods" . 

Two examples of prior art techniques for estimating 
the impulse response of a communication channel are 
10 mentioned. 

The first technique uses particular training 
sequences referred to as constant amplitude zero 
autocorrelation (CAZAC) sequences. These sequences are 
described in an article by A. MILEWSKI : "Periodic 
15 sequences with optimal properties for channel estimation 
and fast start-up equalization", IBM Journal of Research 
and Development, Vol.27, No . 5 , Sept. 83, pages 426-431. 

The GSM cellular mobile radio system uses training 
sequences TS made up of 2 6 symbols ao to 3.^3 taking the 
20 value +1 or -1. These sequences have the following 
properties : 

20 

i=S 
20 

E^i^-k=0 if 0<|k|<5 

i=5 

25 Letting d denote the dispersion of the channel, 

which takes the value 4 in GSM, the estimate of the 
impulse response takes the form of a vector X with five 
components Xq to X4 . 

The received symbol sequence S corresponding to the 

30 training sequence TS is also made up of 2 6 symbols, 

denoted Sq to S25 . The natural assumption is made here 
that the transmitter and the receiver are perfectly 
synchronized, in which case the estimate of the impulse 
response X is given by the following expression: 
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^ 20 
16 X x+k 



for 0<k<4 



The CAZAC technique has the advantage that it is 
very simple to implement. However, it should be noted 
that each component of the impulse response is 
established from only 16 received symbols. Because the 
training sequence is made up of 2 6 symbols and the 
channel dispersion value is 4, there is information in 
the received signal that is not taken into account and 
this degrades performance compared to the theoretical 
ideal . 

The second prior art technique uses the least 
squares criterion. It is described in particular in 
patent applications FR 2 696 604 and EP 0 564 849. It 
uses a measurement matrix A constructed from a training 
sequence TS of length n. The matrix has (n-d) rows and 
(d+1) columns, where d again represents the dispersion of 
the channel. The item in the ith row and the jth column 
is the (d+i-j)th symbol of the training sequence: 



The training sequence is chosen so that the matrix 
a'^A, where the operator represents transposition, 
cannot be inverted. This is inherently the case for 
CAZAC sequences but is also the case for other sequences. 
2 5 The first four symbols Sq through S3 in the sequence 

of received symbols are ignored because they also depend 
on unknown symbols transmitted before the training 
sequence, given that the value of the channel dispersion 
is 4. At the risk of using a misnomer, the received 



4 



signal will therefore be defined as a vector S whose 
components are the received symbols S4, S5, Sg, S25 . 

The estimate of the impulse response then takes the 
following form: 
5 X = (A^A) "V . S 

This least squares technique is slightly more 
complex than the preceding technique but it should be 
noted that the matrix {A^A)'^A^ is calculated only once. 
Note also that each component of the estimate of the 
10 impulse response X is obtained from 22 received symbols, 

rather than 16 as in the CAZAC technique. Improved 
performance can therefore be expected. 

However, regardless of the technique used, the 
impulse response of each channel of the communication 
15 path is considered to be independent of the others. 

A first object of the present invention is therefore 
to provide a method of estimating a communication path 
which takes into account the fact that the various 
antennas are spatially linked. 
2 0 The method in accordance with the invention of 

estimating a communication path is used when the path is 
formed of a plurality of channels and the method 
necessitates an estimate of the impulse responses C^, C2, 
Cn of said channels. The method includes the 
25 following steps: 

- acquiring a space statistic of the transmission path, 

- establishing a corrected impulse response at least by 
weighting said impulse response estimates by means of 
said space statistic and an estimate of the additive 

30 noise of said channels. 

The space statistic advantageously corresponds to an 
estimate of the correlation of the communication channels 
taken two by two . 

In a preferred embodiment of the invention the 
35 estimate of the correlation of the communication channels 
takes the form of a space correlation matrix in which the 
element in the ith row and the jth column is obtained by 
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smoothing the product of the Hermitian transposition of 
the estimated impulse response of the ith channel and the 
estimated impulse response of the j th channel. 

According to an additional feature of the invention 
5 a signal S received by a channel corresponds to a 

transmitted training sequence and the estimate of the 
additive noise N^-^, of that channel is obtained by 
normalizing the energy of the vector (S - AC]^) where A is 
the measurement matrix associated with said training 
10 sequence. 

The normalization can followed by an averaging step. 

Also, if a noise matrix N is formed from the 
estimated additive noise Nq-^, Nqj , of the channels 

and a space -weighting matrix G' is defined on the basis 
15 of said spatial correlation matrix G and said noise 

matrix G' = G(G + N) said corrected impulse responses 
C'l/ C' 2f •••/ C'n are obtained from the following 
expression : 



^ 1 














= G' 




Ic '^V 







The method of estimating a communication path is 
therefore based on the estimate of the impulse response 
of the various channels considered as independent 
channels. Estimation errors are inevitable. As a 

25 general rule, determining the impulse response of a 

single communication channel is a problem that cannot be 
solved exactly in the presence of additive noise. Also, 
the prior art techniques implicitly assume that the 
impulse response can take any form. 

30 Accordingly, a second object of the invention is to 

provide a method of estimating the impulse response of a 
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communication channel which has improved response to 
additive noise, in other words which leads to an error 
lower than the estimation error of prior art techniques. 
This method is advantageously applied to at least one of 
5 the channels forming the communication path before 

establishing impulse responses corrected by weighting by 
means of the space statistic and an estimate of the 
additive noise of the channels. 

According to the invention, this method requires a 
10 signal received by a channel and corresponding to a 

transmitted training sequence. The method includes the 
following steps: 

- acquiring a time statistic of the transmission channel, 

- establishing the estimate of the impulse response of 
15 said channel, which estimate is weighted by said time 

statistic of the channel by means of said received 
signal . 

The time statistic of the channel represents a value 
of the impulse response prior to acquisition of the 

20 received signal. Said weighting introduces the fact that 
the impulse response related to the received signal has a 
value which is probably closer to that prior value than a 
value very far away from it. Thus the estimation error 
is reduced from the statistical point of view. 

25 The time statistic advantageously corresponds to an 

estimate of the covariance of said impulse response. 

A first variant of the method includes the following 
steps : 

- smoothing said impulse response and orthonormalizing by 
30 means of a transformation matrix W to obtain said 

estimate of the covariance which then takes the form of a 
matrix L' , 

- seeking eigenvectors Vi' and eigenvalues A-i' associated 
with that matrix L' , 

35 - estimating the instantaneous impulse response of the 
channel from said received signal and applying that 
transformation matrix W to form a vector X' , so 
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establishing said weighted estimate Xp : 

where Nq is a positive real number representing the 
5 additive noise of said channel. 

The additive noise can be made equal to the smallest 
of said eigenvalues X^' . 

Each eigenvalue of a subset of said eigenvalues X^' 
having a contribution less than a predetermined threshold 
10 can be forced to the value of said additive noise. 

This reduces the complexity commensurately . 
In a second variant of the method the estimate of 
the covariance takes the form of a matrix R and said 
weighted estimate is established as follows: 
15 Xp = (A^A + NoR"^)"W.S 

where A is the measurement matrix associated with said 
training sequence and Nq is a positive real number 
representing the additive noise of said channel. 

The method can include a step of orthonormalizing 
20 said matrix R by means of a transformation matrix W to 

obtain a new matrix R' , the weighted estimate then taking 
the following new form: 

Xp = W^(I + NgR' 'Va'^.S 
where the matrix A' is equal to product of the 
2 5 transformation matrix W and said measurement matrix A. 

The expression (I + NqR' "■"") is advantageously 
calculated by means of the matrix inversion lemma. 

The present invention will emerge in more detail 
from the following description of embodiments of the 
30 invention which is given by way of example only and with 
reference to the accompanying drawings, in which: 

- Figure 1 is a diagram identifying the principal 
steps of an implementation of a method of the invention 
for estimating a communication path, 
35 - Figure 2 shows a first variant of a method of the 



invention for estimating the impulse response of a 
channel , 

- Figure 3 shows a second variant of a method of the 
invention for estimating the impulse response of a 
channel . 

The invention is described as applied to GSM, 
because GSM has the merit of being well-known to the 
skilled person. Thus GSM is described in the interests 
of clarity, but this must not be taken as limiting the 
invention to this system alone. 

Referring to Figure 1, the method of estimating a 
communication path is applied when the path includes at 
least two communication channels, generally n channels. 
Each channel corresponds to a separate antenna. It is 
assumed that an estimate of the respective impulse 
responses C^, , . . . , of each of the channels has been 
arrived at using any of the available techniques. 

The method first acquires a space statistic of the 
communication path. The expression ^^space statistic" 
refers to a set of data reflecting the behavior of the 
path over a predetermined period referred to as the 
correlation period. Because the various antennas are 
fixed, the signals received at the antennas have some 
degree of correlation. The invention aims specifically 
to exploit this fact to improve the quality of the 
estimate of the impulse response of at least one channel. 
For example, this statistic can be obtained by means of a 
spatial correlation matrix G: 



^V(C^Ci) V(C^C2) V(C^CJ 

ViC^C^) V(C^C^) V(C2C^) 



vV(C^C3_) V(C^C2) ViC^Cj , 

where the operator represents the Hermit ian 
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transposition . 

The square matrix G with dimensions (n,n) can 
therefore be represented generically by the element g^j in 
the ith row and the j th column: 
Gij = V(c\Cj) 

The element is obtained by smoothing the product 
c'^iCj using the estimated impulse responses C^, Cj of the 
ith and j th channels obtained during the correlation 
period. This smoothing is an estimate of the correlation 
of the two channels. 

Here "smoothing" is to be understood in a very 
general sense, meaning any operation for smoothing or 
averaging the product c'^iCj over the correlation period. 

A first example of smoothing consists in averaging 
the product over the correlation period: 



The correlation period is assumed to include e 
successive estimates of each of the impulse responses C-^, 
C2, C,. 

A second example of smoothing consists in updating, 
at the pth estimate received for each of the ith and j th 
channels, the smoothing expression Vp _ i(c\Cj) obtained at 
the (p - l)th estimation by means of a multiplier 
coefficient a which has a value from 0 to 1 and is often 
referred to as the smoothing forget factor: 
Vp(c\Cj) = ac\Cj + (l-a)Vp . ^(C^Cj) 

Initialization can be effected by any means, in 
particular using the first estimate obtained or an 
average for the first estimates received obtained as in 
the first example. 

The estimation method then proposes to estimate the 
additive noise Nqi, , present in each of the 

channels by means of the estimates of the respective 
impulse responses C-^, , of the channels. 

Various solutions for obtaining this estimate of the 
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noise are proposed and the case of a single channel, for 
example the first channel, is discussed, given that the 
same solutions apply to each channel . 

A simple solution assigns Nqi a predetermined value 
5 which reflects a threshold below which it is considered 

to be unlikely that the additive noise can descend. This 
value can be determined by measuring a signal to noise 
ratio or from the performance of the receiver, for 
example . 

10 Also, the noise Nq-^ can be estimated by means of the 

estimate of the impulse response of the first channel 
and the corresponding received signal S. Thus the noise 
can be estimated before applying the space-weighted 
estimation method, regardless of which technique is used. 

15 Nevertheless, if this were not the case, there is 
proposed here a method which is appropriate if the 
estimate of the impulse response is acquired using the 
least squares technique. 

It is therefore a matter of estimating the additive 

2 0 noise directly from the received signal S and the 

measurement matrix A. Letting N-^ denote the noise vector 
affecting the received signal : 
S = ACi + Ni 
The estimate of the noise Nqi takes the form: 

25 

Noi=(— ) {S-ACi)'"(S-ACi) 
22 

because the vectors S and N-^ have 22 components . 

This estimate of the additive noise Ngi can 
naturally be averaged or smoothed. 

3 0 A noise matrix N in which n^j denotes the element in 

the ith row and the jth column is constructed from the 
estimates of the additive noise Nqi, N02 , • • '^cm of the 
various channels. 

If i is equal to ±, the element n^j takes the value 
35 of the estimated noise N^^. 

If i and j_ are different, the element n^j is a null 
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element . 

Also, and again in the interests of simplification, 
it is possible to calculate the average value Ng of the 
estimated noise Nqi, Nqj, Nq^ and to force each 

5 element n^j on the diagonal of the matrix N to that 

average value Nq . Letting I denote the identity matrix, 
it follows that the noise matrix N takes the following 
form: 

N = NqI 

10 The correlation matrix G and the noise matrix N are 

used to define a new matrix, namely the space weighting 
matrix G' : 

G' = G(G + N) 

The impulse responses C'^, C'2, C'^ corrected by 

15 means of the space weighting matrix G' are therefore 
defined as follows: 



where . again represents the transposition operator. 
20 At least one of the corrected impulse responses C'l 

is used instead of the estimate of the impulse response 

in the receiver. 

The invention offers the option of improving the 

estimate of the impulse response C^, C^, of each 

2 5 channel before applying the space-weighted method of 

estimating a communication path, i.e. before establishing 

the corrected impulse response (s) C'-^, C'2, C'^- 

A time statistic is therefore acquired for at least 

one of these channels, for example the first channel. 
30 The expression "time statistic" refers to a set of data 

reflecting the behavior of the channel concerned. 













= G' 
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independently of the other channels and over a 
predetermined period referred to as the analysis period. 
It is therefore a representation of the average behavior 
of the channel during the analysis period. This 
5 statistic can be established by any means and anywhere. 

The statistic can be established in any equipment unit of 
the radio communication network. What is important is 
that the receiver is able to acquire this statistic. 

A time statistic of this kind can be obtained in the 
10 following manner, for example. 

Using a method known in the art, an estimate X of 
the impulse response is calculated from the signal S 
received during the analysis period. 

If the least squares technique is adopted, the value 
15 of this estimate X is: 

X = (A^ A) At . S (1) 
It must be remembered that the transmitter and 
receiver are assumed to be synchronized to within better 
than half a symbol, in which case the received signal is 
2 0 the vector S whose components are the received symbols S4 
to S25 synchronous with the symbols a.^^ to a.2s of the 
training sequence TS . Several solutions are available 
for acquiring this synchronization, if not acquired 
already, and two examples of these will be mentioned. 
25 The first solution consists in advancing or delaying 

the received signal by ± symbol periods so that 

Sj^ = (S4_j, S5_j, Sg.j, . . . , S25-j) . 

The estimate Xj is then calculated for each vector 
Sj and the value for which Xj'^.Xj is a maximum is 
30 adopted. This value gives the expected synchronization 
and it is sufficient to replace the vector S in equation 
(1) with the vector Sj„. 

The second solution artificially increases the 
dispersion d of the channel by a predetermined quantity 
35 2q. A modified measurement matrix A^ can then be defined 
with (n-d-2q) rows and (d+2q+l) columns. Assigning n and 
d the respective values 26 and 4: 
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a4+2q....a, 
as +2q....a, 

ag +2, 
a^ +2, 



It is then necessary to reduce the number of 
components of the received signal S by the same quantity 
2q and by convention the modified vector is retained: 



(s 



5+q/ 



- / « 25-q,> 

therefore obtained: 



A modified estimate 
The modified estimate X^, has d+2q+l components: 

X-L , . • ■ / X4 , . . . / 

With the operator . representing the complex 
conjugate, the value of j. between -q and +q which 
maximizes the following expression: 



15 ^x;^,,Xj,k 

k=0 

is then looked for. 

The value j„ determines the estimate X of the 
impulse response for a dispersion d=4 : 

X = (Xjjj, XjM+i, • • • / X-jf^^4) 

20 The synchronization is deduced immediately by 

applying the offset to the received signal S. 
Equation (1) can then be applied again. 
A smoothing matrix L is then constructed by 
smoothing the various estimates X obtained during the 

25 analysis period to obtain an estimate of the covariance 
associated with that impulse response. Here "smoothing" 
is to be understood in a very general sense, meaning any 
operation for smoothing or averaging the impulse response 
over the analysis period. This yields a statistical 
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representation of the behavior of the transmission 
channel. This smoothing can be achieved by either of the 
two methods proposed above, the simplest expression of 
the smoothing matrix L, where m corresponds to the number 
5 of training sequences over which the smoothing is 
calculated, being as follows: 




It is assumed here that the smoothing matrix can be 
10 approximated by the following equation: 
L « (A^A) 'X + R (2) 
in which Ng again represents the noise present in the 
communication channel or additive noise and R is a matrix 
that is usually referred to as the a priori statistical 
15 matrix of the channel because it represents the behavior 
of the channel ignoring noise. 

It is also assumed that the measurement matrix A is 
properly conditioned, i.e. that the eigenvalues of the 
matrix A*^A are very close to each other. In this case, it 
20 is beneficial to normalize the vectors consisting of the 
measurement matrix A, but this must not be seen as 
limiting the invention. 

For this purpose, a transformation matrix W is used 
such that : 
2 5 A' = AW and A'^A' = I 

where I represents the identity matrix. 

Letting L' denote the matrix so defined: 
L = WL'W'', 

it is found that equation (2) can now be written: 
30 L' « NqI + R' (3) 

Note that in a first variant shown in Figure 2 the 
eigenvectors v^' of L' and v^ of R' are identical whereas 
the eigenvalues X^' of L' and X± of R' are offset by Nq . 
Taking the same value of 4 for the dispersion of the 
35 channel, for any i from 0 to 4 : 
Vi' = Vi 
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It is therefore apparent that the eigenvectors and 
eigenvalues of R' and L can be determined in exactly the 
same way, provided that Nq is known. 
5 The step of estimating the noise will be described 

later to clarify the explanation, although it precedes 
the step explained next . 

It is therefore now necessary to seek eigenvalue/ 
eigenvector pairs for the matrix L' or R' . This step 

10 will not be described in detail because it is well-known 
to the skilled person. Moreover, it goes without saying 
that eigenvalues whose contribution is deemed to be 
insignificant can be eliminated. For example, if the 
eigenvalues are listed in decreasing order, the lowest 

15 values whose sum is below a predetermined threshold can 
be eliminated. 

The next step is to estimate the instantaneous 
impulse response X from the received signal corresponding 
to the last training sequence received and using any 

20 technique known in the art. Using the notation X = WX' , 
the latter estimate is weighted by the following method 
to obtain a temporal weighting Xp of the instantaneous 
impulse response: 

To obtain the weighting Xp it is therefore necessary 
to estimate the additive noise Nq . 
30 The noise can be estimated using any of the methods 

referred to above. 

Another possibility is to consider the last 
(smallest) eigenvalue of the smoothing matrix L as equal 
to Nq: 
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= No or A,4 = 0 
Whatever method is adopted, the time weighting Xp of 
the estimate of the instantaneous impulse response can 
then be obtained as indicated above. 

In a second variant, shown in Figure 3, the weighted 
estimate Xp is established directly as: 

Xp = (A^A + NoR"^)"W.S 
Or, using the transformation matrix W defined above: 
Xp = W(I + NqR'"^)"^ A'^.S (4) 
From equation (3) : 

R' = L' = NqI 
Again the additive noise Nq must be estimated. 
An advantageous solution to obtaining the temporal 
weighting Xp is to use the following method. 
The matrix R' is divided by Nq : 



No 

It follows that: 

I + NqR'"^ = I + B'^ 

The matrix inverting lemma is used to calculate the 
weighting matrix P = (I + B'^)'^. 

Accordingly, denoting the canonic vectors e^, the 
following iteration is performed: 

- initialization: 

P = B 

- for ± varying from 0 to d (4 in this instance) : 

Pe.(Pe.)^ 

P = P , 

l + e^Pe^ 

Because P is known, all that remains is to establish 
the weighting Xp from equation (4) . 

Note that the weighting matrix P is not necessarily 
calculated as each new training sequence is transmitted. 
It can be calculated at a slower rate because it varies 
at substantially the same rate as R' and thus more slowly 
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than the received signal S . 

Note also that the weighted estimate is obtained 
without recourse to the instantaneous impulse response. 
It is produced directly from the received signal S. 

Whichever variant is adopted, it is therefore the 
weighted estimate Xp that is advantageously used as the 
estimate of the impulse response to implement the 
space-weighted method of estimating a communication path, 
i.e. to produce one or more corrected impulse responses. 
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CLAIMS 

1. A method of estimating a communication path formed 
of a plurality of channels, the method necessitating an 
estimate of the impulse response C^, C2 , . . ., of said 

5 channels, characterized in that it includes the following 
steps : 

- acquiring a space statistic of the transmission path, 

- establishing a corrected impulse response (C^, C 2, 
. . . , C'n) at least by weighting said impulse response 

10 estimates (C^, C2 , C-J by means of said space 

statistic and an estimate of the additive noise (Nqi, N02/ 
. . . , Nqh) of said channels. 

2. A method according to claim 1, characterized in that 
said space statistic corresponds to an estimate of the 

15 correlation of said communication channels taken two by 
two . 

3. A method according to claim 2, characterized in that 
said estimate of the correlation of the communication 
channels takes the form of a space correlation matrix (G) 

2 0 in which the element (gij) in the ith row and the j th 

column is obtained by smoothing the product {c\Cj) of the 
Hermitian transposition of the estimated impulse response 

(C^) of the ith channel and the estimated impulse response 

(Cj) of the jth channel. 

25 4. A method according to claim 3, characterized in that 
if a signal S received by a channel corresponds to a 
transmitted training sequence the estimate of the 
additive noise (Nq) of that channel is obtained by 
normalizing the energy of the vector (S - AC^) where A is 

3 0 the measurement matrix associated with said training 

sequence . 

5. A method according to claim 4, characterized in that 
said normalization is followed by an averaging step. 
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6. A method according to any of claims 3 to 5 , 
characterized in that if a noise matrix (N) is formed 
from the estimated additive noise (Nd, Nqs, No^) of 

the channels and a space-weighting matrix (G' ) is defined 
5 on the basis of said spatial correlation matrix (G) and 
said noise matrix G' = G(G + N) said corrected impulse 
responses (C'l, C 2^ C'^) are obtained from the 

following expression: 
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= G' 













7. A method according to any of claims 1 to 6, 
characterized in that, if the signal (S) received by a 
channel corresponds to a transmitted training sequence, 

15 the method includes the following steps before 

establishing said corrected impulse response (C'l) of that 
channel : 

- acquiring a time statistic of the transmission channel, 

- establishing the estimate (Xp) of the impulse response 
20 (C-^) of said channel, which estimate is weighted by said 

time statistic of the channel by means of said received 
signal (S) . 

8. A method according to claim 7, characterized in that 
said time statistic corresponds to an estimate of the 

25 covariance of said impulse response. 

9. A method according to claim 8, characterized in that 
it includes the following steps: 

- smoothing said impulse response and orthonormalizing by 
means of a transformation matrix W to obtain said 
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estimate of the covariance which then takes the form of a 
matrix L' , 

- seeking eigenvectors (v^' ) and eigenvalues (X^' ) 
associated with that matrix L' , 
5 - estimating the instantaneous impulse response of the 
channel from said received signal (S) and applying that 
transformation matrix W to form a vector X' , so 
establishing said weighted estimate (Xp) : 

10 Xp =;^|^^^i^(v^x■)jwv^^ 

where Nq is a positive real number representing the 
additive noise of said channel . 

10. A method according to claim 9, characterized in that 
said additive noise (Nq) is made equal to the smallest of 

15 said eigenvalues ) . 

11. A method according to claim 9 or claim 10, 
characterized in that each eigenvalue of a subset of said 
eigenvalues {Xj_' ) having a contribution less than a 
predetermined threshold is forced to the value of said 

20 additive noise (Nq) . 

12. A method according to claim 8, characterized in that 
said estimate of the covariance takes the form of a 
matrix R and said weighted estimate (Xp) is established 
as follows: 

25 Xp = (A^'a + NqR"^) "^A^.S 

where A is the measurement matrix associated with said 
training sequence and Nq is a positive real number 
representing the additive noise of said channel. 

13. A method according to claim 12, characterized in 

30 that it includes a step of orthonormalizing said matrix R 
by means of a transformation matrix W to obtain a new 



matrix R' , the weighted estimate then taking the 
following new form: 

Xp = (I + NqR' "^W^A'^'-S 
where the matrix A' is equal to product of the 
transformation matrix W and said measurement matrix A. 

14. A method according to claim 13, characterized in 
that the expression (I+NgR'"^)" is calculated by means 
the matrix inversion lemma. 
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ABSTRACT 

SPACE -WEIGHTED COMMUNICATION PATH ESTIMATION 

5 The invention relates to a method of estimating a 

communication path formed of a plurality of channels, the 
method necessitating an estimate of the impulse response 
Ci, C2, of said channels, characterized in that it 

includes the following steps: 

10 - acquiring a space statistic of the transmission path, 
- establishing a corrected impulse response (C'-j^, C'2/ 
C'n) at least by weighting said impulse response 
estimates (C^, C2, C^) by means of said space 

statistic and an estimate of the additive noise (Nqi, Nqj/ 

15 . . . , Non) of said channels. 
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